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Dynamical and spectral Dirac systems: 
response function and inverse problems 

A.L. Sakhnovich 


Abstract 

We establish simple connections between response functions of the 
dynamical Dirac systems and A-amplitudes and Weyl functions of the 
spectral Dirac systems. Using these connections we propose a new and 
rigorous procedure to recover a general-type dynamical Dirac system 
from its response function as well as a procedure to construct explicit 
solutions of this problem. 

MSC(2010): 34B20, 35Q41, 35B30, 37D99, 70Q05. 

1 Introduction 

Dynamical systems and corresponding control problems are of great inter¬ 
est and three books on dynamical systems [33, 34, 37] appeared already in 
2015 (see also interesting references therein). In particular, hyperbolic dy¬ 
namical systems (to which class the dynamical Dirac system belongs) are 
important. For dynamical systems and related inverse problems (including 
inverse problems for dynamical Schrodinger and Dirac systems) see, for in¬ 
stance, [6,8,9,12,30,41,46] and numerous references therein. At the same 
time inverse problems for the classical or spectral (frequency domain) self- 
adjoint Dirac system, also called AKNS, ZS or Dirac type system, had been 
actively studied since 1950s (see [38,42]). Various interesting results were 
published last years (see, e.g., [1,2,16,20,22,35,45,47,54]). 
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The interconnections between dynamical and spectral Dirac systems, 
which we establish here, provide a new and rigorous way to recover a general- 
type dynamical Dirac system from the response function and to construct 
explicit solutions of this problem as well. These interconnections open also 
further possibilities for the study of both dynamical and spectral Dirac sys¬ 
tems. 

Spectral Dirac system has the form 

y\x,z) = \{zi + iV{x))y{x,z), a; > 0 (^y\x, z) := ^{x, z)^ , (1.1) 

where 

, mi -|- 1712 =■ (1-2) 

Here is the m^ x m* identity matrix and v{x) is an mi x m 2 matrix 
function. 

As we already mentioned, dynamical Schrodinger and Dirac systems as 
well as their connections with response functions and boundary control are 
of growing interest [5,6, 8, 9]. Dynamical Dirac system (Dirac system in 
the time-domain setup) was studied in the important recent paper [9]. The 
dynamical Dirac system considered in [9] is an evolution system of hyperbolic 
type and has the following form: 


J = 


‘■mi 

0 


‘■m2 


D = 


iut Jux + Vu = 0 (a; > 0, t > 0); 


u = 

'ui 

, J = 

1 1 
1 0 

, V = 

1 - 


“2. 



Iq -pj 



( 1 . 3 ) 

( 1 . 4 ) 


where p = p{x) and q = q{x) are real-valued functions of x, and initial¬ 
boundary conditions are given by the equalities 


M(a;,0) = 0, a; > 0; Mi(0,t) = fit), t>0. (1.5) 


Here / is a complex-valued function (so called boundary control) and the 
input-output map (response operator) R : Mi(0, •) —)■ ^ 2 ( 0 , •) is of the con¬ 
volution form Rf = if r * f. The inverse problem consists in recovery 
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of the potential V from the response funetion r. This inverse problem was 
considered in [9] using boundary control methods. 

We note that recent results on dynamical Schrodinger and Dirac equa¬ 
tions are based on several earlier works. In his paper [10] from 1971, A.S. Bla- 
govescenskii considered dynamical system 

titt T d (^*d) 

with boundary control u{0,t) = f{t), and solved inverse problem to recover 
Q from /. A.S. Blagovescenskii established important connections between 
his problem and spectral theory of string equations. This work was developed 
further in [7] (see also references therein), where response operator appears in 
inverse problem. Finally, the inverse problem to recover the matrix potential 
Q{x) of the dynamical Schrodinger equation 


utt - Uxx + Q{x)u = 0 (1.7) 

from the response function was considered in [4] by S. Avdonin, M. Belishev, 
and S. Ivanov. 

Similar to the case of the spectral Dirac and Schrodinger equations, 
the dynamical Dirac equation is a more general object than the dynamical 
Schrodinger equation. More precisely, setting in (1.3) 

p{x) = 0, q{x) = gx{x)/g{x), where gxx{x) = Q{x)g{x), (1.8) 

and rewriting (1.3), (1.4) in the form 

{ui)t = i{{u 2 )x + qu^) , {u 2 )t = i( - {ui)x + qui) , (1.9) 

we obtain a dynamical Schrodinger equation 

(^i)tt (^i)xx i^qx "F q )^i (^i)xx (^9xx/g^ui 

= {Ul)xx-QUl- ( 1 - 10 ) 

For interesting applications of the interconnections between spectral Dirac 
and Schrodinger equations see, for instance, the papers [11,20,21,26,27] and 
references therein. 
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In our paper, we come from the dynamical to the spectral Dirac system 
taking Fourier transformation of both parts of (1.3). For that we should be 
able to estimate the behavior of solutions u of (1.3)-(1.5). Our next section 
is dedicated to the necessary estimates on u. A simple connection between 
response function of the dynamical Dirac system and Weyl function of the 
corresponding spectral Dirac system is derived in Section 3. Recovery of the 
dynamical Dirac system from the response function is described in Corollary 
4.5 in Section 4. The connection between response function and so-called 
A-amplitude is described in Section 4 as well. We note that interrelations 
between dynamical and spectral Dirac systems provide a procedure for solv¬ 
ing the general-type inverse problem for dynamical Dirac system in much 
greater generality than in [9]. 

Explicit methods in spectral theory and in construction of solutions of 
equations is an interesting domain, which is actively developed (to a great 
extent independently from the general theory). In particular, the methods 
include different versions of Backlund-Darboux transformations, Crum-Krein 
transformations and commutation methods (see, e.g., [14,15,17,19,25,29,36, 
40,44,51,56] and various references therein). In Section 5, we recover ex¬ 
plicitly dynamical Dirac system from response function using GBDT (gener¬ 
alized Backlund-Darboux transformation) technique from [32,49,51,54] (see 
also references therein). 

In order to avoid the usage of the same letters for different items, we often 
use the calligraphic font (e.g., V) for functions and operators corresponding 
to dynamical system (1.3). As usual, M stands for the real axis, C stands 
for the complex plane, C+ is the open upper half-plane {z : A (z) > 0} 
and Cm is the half-plane {z : A( 2 :) > M > 0}. The notation 6 stands for 
the complex conjugate of 9 when 0 is a complex number, and for the vector 
with the entries which are complex conjugates of the entries of 9 when 9 is 
a vector. By l.i.m. we denote the entrywise limit of matrix functions in the 
norm on hnite intervals, semiaxes or R depending on the context. We say 
that the matrix function is boundedly (continuously) differentiable when its 
derivative is bounded in the matrix norm (continuous). By B{H) we denote 
the class of bounded linear operators acting in the Hilbert space H, I is the 
identity operator, Im is the m x m identity matrix and A* stands for the 
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Hermitian adjoint of an operator (or a matrix) A. 


2 Preliminaries and estimates on the solu¬ 
tions u of the dynamical Dirac systems 

According to [9, Theorem 1], in the case where p, q, f are continuously differ¬ 
entiable (i.e., p, q, f E C^) and /(O) = /'(O) = 0, there is a unique classical 
solution M of (1.3), (1.5) and this solution admits representation 


u = u{ A w-^; 

ui{x,t) = f{t - x) J , /(t) = 0 (t<0); 

(2.1) 

w^{x, t) = w{ 

[x,t) = [ f{t — s)K{x,s)ds {t>x>0), 

J X 

(2.2) 

w{x, t) = 0 

{x > t > 0), 

(2.3) 


where k{x, s) [x < s) is continuously differentiable. In particular, formulas 
(2.1) and (2.3) yield: 

u{x,t) = 0 for 0 <t < X (hnitenessof the domain of influence). (2.4) 

Representation (2.1)-(2.3) is proved in [9] using Duhamel formula. Moreover, 
the second term w-^ = w in the representation (2.1) of u admits (see [9, p. 6]) 
an expansion 

w{x,t) = for t>x>0, (2.5) 

A:>0 

{Ag) (x, t) = -S{V{x)g{x, t )), (2.6) 
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and the operator S is given by the right-hand side integrals in [9, formulas 
(1.8)-(1.11)]. Namely, for t > x we have 


{Sh)^{x,t) 


{Sh)^{x,t) 


r*A 2 


2^/2 Oai 


fAs 


(ihi -|- h2)M — / (i/ll + h2)M 


' A 2 


r*A 4 


-h / (i/ll - h 2 )di ) , h{x,t) = 
J\i 


hi{x,t) 
[h2{x,t)\ ’ 


2^/2 


r*A 2 


{hi — ih2)di 


fAs 


{hi — ih2)di 


' Xi 


' A 2 


.A 4 


(/ii + '\h2)M J , 


'Ai 


(2,7) 


( 2 , 8 ) 


where Aj G for i = 1, 2, 3,4 and Ai = (x, t), X 2 = (0, t — x), A 3 = (/ — x, 0), 
A 4 = (x-l-/, 0). Here -di stands for the integral along the interval [A*, A^] 
and i is the length. 

If h{x,t) = 0 for t < X, we have [Sh){x,t) = 0 for t < x, and therefore, 
in view of ( 2 . 1 ), we shall need here only the expression for Sh when t > x. 
Moreover, since we assume that /(/) = 0 for / < 0, it follows from (2.3) and 
(2.6) that (2.5) holds for / < a; as well. 

Let us also assume that V, / and /' are bounded: 


sup ||V(a;)|| < Ml, 

a ;>0 


sup 

l>0 



< Co, 


sup 

l>0 


fit) 


1 

i 


< Co. 


(2.9) 


Proposition 2.1 Let p, q, f be continuously differentiable and let equalities 
f{t) = f'{f) = 0 hold fort < 0. Assume that (2.9) is valid. Then the solution 
u of the dynamical Dirac system (1.3), such that (1.5) and (2.4) are valid, 
satisfies the following inequalities 


||M(a;,/)|| < coe^*, ||tit(a;,/)|| < coe^*, \\uXx,t)\\ < M2e^\ (2.10) 

where a; > 0 and t > 0, M 2 > 0 is some constant, and M = 2\/2Mi. 

Moreover, the functions u{x,t), Ut{x,t) and Ux{x,t) are continuous in 
the quarterplane a; > 0, / > 0. 


Proof. Using (2.6)- (2.8) and the hrst inequality in (2.9), it is easy to show 
that 


\\{Ag){x,t)\\ <2V2ckMit’^+^/{k + l) 
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for the case that || 5 f(a;, t)|| < t>x (recall that we are interested in the 
case g{x,t) = 0 for t < x). Hence, relation (2.1) and the second inequality 
in (2.9) imply that 

||(^"+V)(a;,t)|| <co(2 V2Mi)'^+V+V(A; + l)!. (2.11) 

It follows from (2.5) and (2.11) that 

lk(a;,t)|| < Co(e^* - 1) {M = 2y/2Mi). (2.12) 

Formulas (2.1), (2.9) and (2.12) yield the inequality ||M(a;,t)|| < Coe^h 

It is easy to see that the functions {x, t) {k > 0) are continuous 

at X = t and = 0. Thus, using the hrst and third inequalities in 

(2.9) we (similar to the derivation of formula (2.11)) derive 

||(.4'=+‘u{),(x,t)|| <5j(2V2M,)‘+'t'=+V(l-+l)!, (2.13) 

and the second inequality in (2.10) follows. Finally, the third inequality in 

(2.10) is immediate from (1.3), the hrst inequality in (2.9) and the hrst and 
second inequalities in (2.10). 

Recall that u{x,t) is continuous at x = t. According to (2.1) and (2.2), 
the functions Ux{x,t) and Ut{x,t) are continuous at x = t. It is evident that 
u, Ux and Ut are continuous outside x = t as well. Thus, u, Ux and Ut are 
continuous in the quarterplane. ■ 

Proposition 2.1 yields the following corollary. 

Corollary 2.2 We can apply to u the transformation: 

e^^^u{x,t)dt, z G Cm = {z : ^( 2 ;) > M}, (2-14) 

where u stands for the Fourier transformation of u {and Fourier transforma¬ 
tion is taken, for the sake of convenience, for the fixed values x, z ). Moreover, 
the same transformation can he applied to Ut and we have 

poo 

i / F^^Ut{x,t)dt = zu{x, z), zeCm- (2T5) 

Jo 
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Using the mean value theorem and the fact that is continuous and satishes 
the third inequality in ( 2 . 10 ), we obtain our next corollary. 

Corollary 2.3 Fourier transformation can be applied to and 

F^^Ux{x,t)dt = -^u{x,z) = u'{x,z), z e Cm- (2.16) 

(JjJb 

Now, applying the Fourier transformation to the dynamical Dirac system 
(1.3), we derive 



zu{x, z) + Ju'{x, z) + V(x)u(x, z) = 0. (2.17) 

We note that various generalized Fourier transformations are successfully 
used (see, e.g., [3,18,54]) for solving inverse problems. 

An estimate for k{x, t) is also necessary for our further considerations. 
According to [9, p. 7], the formula 

t) = ^A^ {Aui ); {Aui) (x, t) = - ^lim^ 5 ^(p + ig )4 

5e{t) ;= 1 /e for 0 < t < e, 5e{t) := 0 for t > e 

is valid. Using (2.7)-(2.9) and the second equality in (2.18), it is easy to 
show that 



\\{Aui){xA)\\<‘2V2Mr. (2.19) 

From the hrst equality in (2.18) and from (2.19), similar to the proof of 
Proposition 2.1, we obtain the following proposition. 

Proposition 2.4 Let the conditions of Proposition 2.1 hold. Then 

\\K{x,t)\\ < (2.20) 


3 Response and Weyl functions 

Recall that the m 2 x mi Weyl matrix function (Weyl function) (p{z) of the 
spectral Dirac system (1.1) with the locally summable potential V is uniquely 



determined by the inequality 


\lmi (p{z)*].Y{x,z)*Y{x,z) dx < oo, zeC+, (3.1) 

[}p[z)\ 

where Y(x,z) is the fundamental m x m solution of (1.1) normalized by 
the condition Y{0,z) = Im (see [22], [54, Section 2.2] and some references 
therein). Here C+ is the open upper half-plane. Weyl functions cp introduced 
by (3.1) are analytic and contractive in C+. 

In view of (2.1) and (2.14) we see that 

u{x, z) = f{z)e^^^ ^ +w{x,z), (3.2) 



where (according to ( 2 . 12 )) 

Mx,z)\\= j\^^^w{x,t)dt < (3.3) 

Hence, for 2 : G Cm, the relation w{x,z) G L 2 ( 0 ,cxd) holds and implies 
u{x,z) G L^O, cxd). In other words u is the Weyl solution of the Dirac 
equation (2.17). The transformation of the spectral (i.e., frequency-domain) 
Dirac system (2.17) (and its solutions) into the equivalent form (1.1) is given 
by the formula 

1 ” i 1 ” 

y = }Cu, v = iq-p, /C = ^ ^ , (3.4) 

which means that Ku is the Weyl solution of (1.1), (1-2), where mi = m 2 = 1. 
Finally, in view of (1.5), (2.1) and (2.2) we obtain 

wi(0,t) =/(t), U2{0,t) = if{t) + [ r{t - s)f{s)ds, (3.5) 

Jo 

where r(t) = K,2{0,t) is the response function [9]. Taking into account (3.5) 
and estimates (2.9) and (2.20), we have 

ui{0, z) = f{z), U2{0,z) = f{z){r{z)+i). (3.6) 


9 



Now, using (3.6), we express Km in terms of the normalized fundamental 
soluton Y: 


y{x, z) = Ku{x, z) = f{z)Y{x, z)K 


riz) +1 


(3.7) 


Recall that y = Ku is the Weyl solution, that is, y G 1^2(0 )that 
the Weyl function if[z) is uniquely determined (see dehnition (3.1)) by the 
condition 


Y {x, z) 


1 

_^{z)_ 


e Ll{0,oo). 


Therefore, formula (3.7) implies that 


(p(z) = r{z)/(r{z) + 2 i). 


(3.8) 


Proposition 3.1 The response function r{t) of the dynamical Dirac system 
(1.3) is connected with the Weyl function (p(z) of the corresponding spectral 
Dirac system (1.1) {where mi = m 2 = 1) via equality (3.8). 

We note that, for the classical case mi = m 2 = k, it is often more convenient 
to introduce Weyl functions (Ph{z) (see, e.g., Dehnition 1.51 [54]) which be¬ 
long to Herglotz class instead of being contractive. (Recall that Herglotz ma¬ 
trix functions in C+ are such matrix functions '4){z) that i{f){z)* —ip{z)) > 0 
for all G C+.) Namely, Weyl matrix function (Ph{z) is uniquely determined 
by the inequality: 



i^h{z)* ] KY{x,zyY{x,z)K* 


1 Ik —Ik 

\/2.4 Ik 


z G C+. 


Ik 

-Wh{z) 


dx < 00 , 


(3.9) 


(3.10) 


A simple connection between ip and pn, which easily follows from (3.1) and 
(3.9), implies (for k = 1) that (3.8) may be rewritten in the form 


Ph{z) =r{z)+i. (3.11) 

Thus, our procedure [50, 52, 54] to recover potentials V of Dirac systems 
(1.1) can be used to recover potentials V of dynamical Dirac systems from 
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response functions. Moreover, even for the case of continuously differentiable 
potentials, equality (3.11) generates important new results, especially on 
explicit recovery of potentials. 

Open Problem. Show that equality (3.11) holds for much wider classes 
of potentials V then continuously differentiable potentials and is also valid 
for the non-scalar case mi = m 2 > 1. 

4 Response function, A-amplitude 
and general-type inverse problem 

4.1 Inverse problem for the spectral Dirac system 

Since V is bounded, the procedure to recover locally bounded potentials 
from Weyl functions, which is given in [50, Theorem 5.4], suffices for our 
purposes. A quick summary of this procedure (with some functions and op¬ 
erators multiplied, for convenience, by corresponding constant scalar factors) 
is presented below (and is close to the summary in [24, Section 2]). 

First, introduce a family of convolution operators: 

d C’' 

S'/ = — / s{x — t) ■ dt, s{x) = —s{—x)*, (4.1) 

dx Jq 

where the index in the notation of the operator S'/ indicates the space 
L^(0,/), in which S'/ is acting, and the matrix function s{x) is associated 
with Dirac system. More precisely, s is introduced (for a; > 0) via the Weyl 
function ipn- 

^ (^^e^"^l.i.m.a^oo j + ir/)“Vr/(^ + ■ (4-2) 

Here p > 0 and (pn is the Weyl function of a spectral Dirac system (1.1), 
where mi = m 2 = k and V is locally bounded on [0, cxd). Recall that ipn is 
uniquely determined via the inequality (3.9). The notation l.i.m. denotes in 
(4.2) the entrywise limit in the norm of T^(0, cxd). It was proved in [50] that 

a + ih)"V(^ + ih) e LLfc(-oo, 00) 
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for every fixed t] > 0, that l.i.m. on the right-hand side of (4.2) is differen¬ 
tiable (and does not depend on 77 > 0), and so s{x) is well-defined. Moreover, 
s{x) is boundedly differentiable on the intervals ( 0 , /), <s(-|- 0 ) = 4/^, and the 
operators Si are bounded and positive (i.e., Si > 0) for all / > 0. Thus, we 
have 

Si = I + f uj{x — t) ■ dt > 0 , u{x) := s'{x), (4.3) 

Jo 

uix) = ijji—x)*^ s' := (4.4) 

dx 

and Si, Sf^ e B(^Ll{0,l)y 

Next, introduce the matrix-functions 

e,{x) = [h 0]Y{x,0)K\ 02 (a:) = [O h]Y{x,0)K*. (4.5) 

According to [50, formula (4.16)] we have 

^ 2 ( 3 :) = ^ ([-4 h]d?j , (4.6) 

where Sf^ (/ = 2x) is applied to [2s(t) 4] columnwise. Recall that Y 

satishes (1.1) and is normalized by Y{0,z) = Im, and that K is given by 
(3.10). Hence, 

KY{x,0)*jY{x,0)K* = J, Y{x,0)K*JKY{x,0)* = j, J := ^ ■ 

[4 0 

(4.7) 

Using (1.1), (4.5) and the second relation in (4.7), we obtain 

v{x) = i 6 '^{x)J 62 {x)*. (4.8) 


Finally, from (1.1), (3.10), (4.5) and the second relation in (4.7) we derive 
the equalities 


0i(O) 



9i{x)J92{x)* = 0 , 


9'^{x)J9i{x)* = 0, 


(4.9) 


which uniquely determine 9i, assuming that ^2 is already given. We can now 
formulate Theorem 5.4 from [50]. 
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Theorem 4.1 Let ipn be the Weyl function of a spectral Dirac system (1.1), 
where mi = m 2 = k and V is locally bounded, that is, 

sup ||l^(a;)|| < cxD for any / > 0. (4-10) 

0<x<l 

Then V can be uniquely recovered from (fn via the formulas (1.2), (4.8) and 
(4.6), (4.9). Here u on the right-hand side of (4.6) is obtained from (4.2) 
{after using u = s') and S 2 X on the right-hand side of (4.6) is given {using 
uj) in (4.3). 

Remark 4.2 For the case mi = m 2 = k = 1, equalities (4.9) are equivalent 
to a much simpler relation 


01 = -92J, (4.11) 

and we recover 9i from (4.11) instead of using (4.9) {if 92 is given). 

In the seminal paper [28] (see also [55] as well as some references in 
[28,55]), the high energy asymptotics of the Weyl functions of Schrodinger 
operators is expressed in terms of the so called A-amplitndes. High energy 
asymptotics of the Weyl fnnctions of self-adjoint (spectral) Dirac systems 
was stndied in [16,50] (see also [48, formula (32)]). According to [50, formula 
(3.25)], the asymptotic equalities 

Th{z) = ilk2i [ [s'{x))*dx o , | 2 :| ^ cxd (/ > 0) (4.12) 

Jo 

hold in all the angles cA(; 2 ) > |3f^(;2)| in C_|_. Formula (4.12) is an analog of 
the asymptotic expression for the Weyl fnnction of a Schrodinger operator in 
terms of the A-amplitnde. 

Remark 4.3 Since 2i (s'(a;))* is an analogue of the A-amplitude for the 
Schrodinger operator case, we call this expression the A-amplitude of the 
spectral Dirac system. We note that, in M. G. Krein’s terminology, u = s' is 
called the accelerant [39] {see further explanations in [2f). 
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4.2 Response function and 74-aniplitude 

Recall that dynamical Dirac system (1.3) generates corresponding spectral 
Dirac system via formulas 


V = 


P Q 
q -p_ 


R = 


0 V 

V 0 


, V = iq — p. 


(4.13) 


Theorem 4.4 Letr{t) he the response function of a dynamical Dirac system 
(1.3) and let the conditions of Proposition 2.1 hold. Then s given by formula 
(4.2) is well-defined and differentiable, and we have the equality 


r{f) = 2is'{f), 


(4.14) 


that is, the response function r of the dynamical Dirac system coincides with 
the A-amplitude of the corresponding spectral Dirac system. 


Proof. We rewrite (3.11) in the form 


Th{z) 


1 = 


P^*r{t)dt, z e 


-M- 


(4.15) 


In view of the equality r(t) = K,2{0,t) and inequality (2.20), formula (4.15) 
yields 

— ^ = ^d.m.a^oo [ D^^'^^^^^r{t)dt, where ?]> Mis fixed (4.16) 
Jo 

and ^ Ihs-l is, hi.m. is taken in the norm L‘^{—oo, oo). Hence, putting 

r{t) = 0 for t < 0, (4.17) 

and taking inverse Fourier transformation, we derive 

(4.18) 


e = l.i.m.„„o^ / e + W) - 


for each hxed p > M. Thus, it is immediate that on each semiaxis {—oo,l) 
we have 


r{t) = l.i.m.a^oo—e"'* / e + ip) - (4.19) 
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By standard calculations, using (4.19), we obtain 




r(t)dtdy 


-lim 

2,71 a —>-00 


-i(€+i„)a. + CiX + ^2 

+ W 


ZTT 


(4.20) 

-iS, + + Cix + C2, 

(^ + w 


where l.i.m. is taken on the intervals (0,1 ) with respect to x. 

On the other hand. Residue Theorem implies that, uniformly for a; > 0, 
the equality 


lim 


3-k^(^ + iri)-^d^ = -2nxe-^^ 


(4.21) 


holds. Hence, formula (4.2) can be rewritten in the form 






(4.22) 


Finally, (4.14) follows from (4.20) and (4.22). 


Corollary 4.5 Letr{t) he the response funetion of a dynamical Dirac system 
(1.3) and let the conditions of Proposition 2.1 hold. Then the function v{x) 
is uniquely recovered from r{t) using the formula 


s{x) = 


1 

2 



a; > 0, 


(4.23) 


and the procedure given in Theorem f.l {more precisely, via the formulas 
(4.8) and (4.3), (4.6), (4.11)). Next, the potential V{x) is uniquely recovered 
fromv{x) using (4.13) or, equivalently, using the first equality in (4.13) and 
the formulas 


p = —3fJ(n), q = 9=(n). 


(4.24) 
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4.3 Extensions and conditions on the 74-aniplitude 
and on response function 

From (4.3) and (4.6) it is clear that in order to recover the values of the 
potential V (of the spectral Dirac system) on [0,/] we need to know u(t) 
or, equivalently, s(t) on the interval [0,2/]. In the paper [9], the dynamical 
Dirac system is considered on the square {{x,t) : 0<x<T,0<t<T}. In 
a similar way, this yields the necessity to introduce an ’’extended” problem 
(and to extend boundary condition / and response function r on the interval 
[0, 2T]) in order to recover V{x) on [0,T] from r. In other words, from r 
given on the interval [0,T], the potential V is recovered on [0,T/2] only. 

It is also interesting to compare the conditions on the Weyl function 
and 4-amplitude for the spectral Dirac system with the conditions on the 
response function for the dynamical Dirac system. For the case that v is 
a square matrix function (or a scalar function), sufficient conditions for a 
Herglotz function ipu iohe a. Weyl function can be given in terms of the 
spectral function [42,43], and the spectral function is connected with ipn via 
Herglotz representation. On the other hand, positive operators Si are also 
recovered from the spectral function (see [54, Theorem 2.11]), and in this 
way the conditions from [42, 43] are related to the condition of positivity 
of Si- The invertibility of the convolution operators (which up to constant 
factors coincide with our operators Si) is required in [2,38] and provides the 
positivity of Si as well. Finally, for the case of the mi xm 2 matrix functions v, 
sufficient conditions for (p{z) to be a Weyl function are local boundedness of 
the 4-amphtude and invertibility (positivity) of the corresponding operators 
Si, see [54, Theorem 2.54]. We see that the positivity of the operators Si 
acting in L^(0,/) is required from (pn to be a Weyl function of a spectral 
Dirac system. The condition for r to be a response function, which is given 
in [9], is again the positivity of a structured operator S acting in the vector 
space L 2 ( 0 , 2 T), although certain discrepancies in the dehnition of S exist 
(see [9, pp. 18 and 23]). 

We note that necessary and sufficient conditions for the solvability of the 
inverse problem, which was studied in [10] for system (1.6), are again formu¬ 
lated in terms of positivity of certain convolution operators (Krein integral 
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operators). 


5 A special type of inverse problem: 
explicit recovery of the potential 

Explicit construction of the Weyl function (direct problem) and explicit re¬ 
covery of the potential V of the spectral Dirac system (1.1) from the Weyl 
function (inverse problem) were dealt with in [23,31,32,51]. We shall use the 
formulation of these results (for the case mi = m 2 = k) from [51, Subsection 
5.1.1]. 

The class of potentials that we recover is called pseudo-exponential [31]. 
Each potential from this class is determined by a hxed number n G N, by an 
n X n matrix A and by two n x k matrices 'di and ^2 such that 

A-A* = iA(0)jA(0)*, A(0) := [di ^ 2 ] . (5.1) 

Then V corresponding to the triple of parameter matrices {A, di, ^ 2 } is 
determined by the second equality in (1.2) and by the formula 

v{x) = -2idte'"^‘^(a;)-^e'"^d2, (5.2) 

where S{x) is given by the equations 

S\x) = K{x)A{x)*, S'(0) = 4; A'(a;) =-iAA(a;)j, A(0) = [di ^ 2 ] . 

The formula dehning S and A above is easily rewritten in the form 

S{x) = In + f A{t)A{t)*dt > 0, (5.3) 

Jo 

A{x) = [Ai(a;) A 2 (a;)] = [e“'®^-di e^*^d 2 ] . (5.4) 

Weyl functions (fn of Dirac systems with pseudo-exponential potentials are 
rational functions belonging to Herglotz class (i.e., \{ipH{.z)* — (Ph{.z)) > 0 
for 2 ; G C+). More precisely, the next statement is valid. 

Proposition 5.1 Let v he a pseudo-exponential potential {i.e., let v admit 
representation (5.2)-(5.4), where (5.1) holds). Then the Weyl function ipn 
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of the spectral Dirac system (1.1), (1.2) with this v (mi = m 2 = k) is given 
by the equality 

(Ph{z) = + a := A - + ^ 92 )*. (5.5) 

Vice versa, each proper rational matrix function (j){z) such that 

lim (l){z) = Ik, i{(l){z)* — (fiz)) >0 for z E C+ (5.6) 

is the Weyl function of a unique spectral Dirac system with a pseudo-exponential 
potential. 

In order to deduce which class of response functions corresponds to Weyl 
functions of the form (5.5) {k = 1), we take into account the equality 

f — C)~^dz = exp{—ixC), (5.7) 

27ri Jr 

which holds [53, p. 557] for anti-clockwise oriented contours F and matrices 
C such that the spectrum of C is situated inside F. In view of (5.5), we 
rewrite (4.19) in the form 

r{t) = {-2i)\A.m.a^ooPz [ e~'''\zln - a)~^dzili, (5.8) 

27ri Jr, 

where Fq are anti-clockwise oriented contours consisting of the points 

{z : jz — irj\ = a, 9^(2:) < ri}U {z : —a < z — if] < a}, 

p > 0 is sufficiently large and hxed, di and ^2 are column vectors (■di,'d 2 E 
C""), and l.i.m. may be considered on any hnite interval (0, /). It is immediate 
from (5.7) and (5.8) that 

r{t) = -2id;e"^‘“di. (5.9) 

We easily check also directly that for r of the form (5.9) we, indeed, have 

r{z) = 2Dl{zIn - a)~^Di (^>(^) > ||q;|1), (5.10) 

that is, 

(Ph{z) =r{z)+i = i + 2dl{zln- (^ 5 ( 2 :) > ||a||). (5.11) 

Thus, Proposition 5.1, relations (3.11) and (5.10) and uniqueness in the re¬ 
covery of V from ifH (see Theorem 4.1) imply the following theorem. 
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Theorem 5.2 Letr{t) he the response funetion of a dynamieal Dirae system 
satisfying conditions of Proposition 2.1 and assume that r(t) admits repre¬ 
sentation (5.9), where the matrix A = a + 101 + -d^)* and column vectors 

dti E {i = 1,2) satisfy (5.1). Then the potential V of this dynamical 
Dirac system is given by the first eguality in (4.13) and formulas (4.24) and 
(5.2)-(5.4). 

We note that identity (5.1) may be rewritten in an equivalent form 

a - a* = -i{'di + id2) + A)*. 

Remark 5.3 The reguirement (in Theorem 5.2) that r should admit rep¬ 
resentation (5.9) could be substituted by some reguirements on the Fourier 
transform r ofr {more precisely, by the conditions (5.6) on (f){z) := r(z) +i). 
In that case 0 is the Weyl function ipn of some spectral Dirac system with a 
pseudo-exponential potential, and so matrices a and Di are recovered from (fn 
following the procedure from [51, Theorem 5.4] (see also references therein). 

It would be of interest also to construct explicit solutions u of dynamical 
Dirac systems similar to the construction of explicit fundamental solutions 
of the spectral Dirac systems in [23,31,32,51]. 
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